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ABSTRACT 

We analyse the motion of single optically thick clouds in the potential of a central mass 
under the influence of an anisotropic radiation field ~ | cos(0)|, a model applicable to 
the inner region of active galactic nuclei. Resulting orbits are analytically soluble 
for constant cloud column densities. All stable orbits are closed, although they have 
non-trivial shapes. Furthermore, there exists a stability criterion in the form of a 
critical inclination, which depends on the luminosity of the central source and the 
column density of the cloud. 



1 INTRODUCTION 

Most - if not all - galaxies harbour supermassive black holes 
(SMBH) at their centres. Gas infall from larger scales leads 
to the formation of an accretion disc, where gravitational 
energy is dissipated to viscously heat the gas. The 
emerging radiation - peaking in the ultraviolet (UV) /optical 
wavelength regime - illuminates the so-called Broad Line 
Region (BLR) on small scales as well as a large gas and 
dust reservoir on parsec scale distance from the SMBH. 
The latter was postulated in order to unify two classes 
of observed objects by means of an inclination effect, 
allowing only unobscured views of the central region for 
face-on orientation of this so-called dusty torus. This is 
the essence of the widely acc epted Unif i ed S cheme of 
Active Galactic Nuclei (AGN, lAntonuccil 1 19931 ). Recent 
observations and modelling efforts provide evidence that 
the inte rnal structure of those t ori is likely c l umpy 
( e - fi- 



nal structure or tnose t ori is likely c lumpy 
Krolik fc Begelmanl Il988l; INenkova et all 120021; 



Honig et all l200d : iTristram et all 120071 ; ISchartmann et all 
20081 ). Successively improved radiative transfer simulations 



have been used to draw conclusions about the properties 
of indivi dual clouds and the geometry of their distribution 
(e. g. INenkova et all 120081 ; iHonig fc Kishimotol l2010l ; 
iRamos Almeida et alj 20091 ). Although differing in detail, 
various models are able to reproduce several types of 
observed spectral energy distributions and interferometric 
data successfully. However, the cloud distribution in these 
models is static and there is only little known about the 
dynamical properties of such cloud systems. Concerning 
the structure of the BLR, several observed occultation 
events suggest the existe nce of a clumpy structure as well 
(e. g. iRisaliti et aTll201ll ). 

In both situations, radiative forces on either gas or dust 
clouds can be of similar strength as gravitational forces. 
The particular significance of anisotropic radiation pressure 
for the distribution of dusty gas has been pointed out by 
ILiu fc Zhangl l|201ll ). In order to understand the dynamical 
stability properties of cloud orbits under the combined 
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effects of strong gravity and radiation pressure, we start 
our analysis by investigating the dynamics of simplified 
single-cloud models. Many aspects of orbits in r adiation 
fields have already been explor ed, for example by ISaslawl 
l|l97Sl) or iMioc fc Radul l| 19921 ). who use a perturbative 
ansatz to examine orbits in general time-varying and 
anisotropic radiation fields, respectively. We aim to analyse 
particularly the orbits of gas and dust clouds in the central 
vicinity of AGN, simplified to a central mass and the 
constant rad iation field of a surro unding accretion disc. In a 
similar setup lKrause et all (|201lf ) investigated perturbations 
to circular orbits of pressure-confined BLR clouds. We 
derive the orbit equation analytically under the alternative 
assumption of constant cloud column density. This will 
enable us to constrain cloud parameters as well as the 
distribution of clouds in AGN tori or the BLR. 

In Section [2] we calculate the effective force acting on 
a cloud and solve the equation of motion in Section [3] We 
proceed to analyse orbit stability in Section [4] and describe 
the characteristics of the orbit family in Section [5] which 
is followed by a critical discussion (Section [5} and the 
conclusions (Section 0. 

2 FORCES ON OPTICALLY THICK CLOUDS 

For the work presented in this Letter, we concentrate on 
the two dominant forces close to a black hole plus accretion 
disc system: the force F gra v on a cloud of mass m in the 
gravitational potential of a central point mass M at distance 
r from the cloud is given by Newton's law: 

GMm 

* grav — 2 ^ r 1 \ ) 

where G is the gravitational constant and e r is the unit 
vector in radial direction. Often of similar dynamical 
importance is the radiation pressure exerted from a viscously 
heated accretion disc. Compared to the cloud distances of 
interest here, the accretion disc emits the bulk of its energy 
close to its inner edge and can hence be considered point-like. 
Assuming that the cloud is optically thick for the whole 
wavelength range where most of the energy is radiated, the 
additional force by radiation pressure can be written as 
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C 

with the anisotropic radiation flux 
L[ cos (0)| 



2-KT 2 



(2) 



(3) 



where a is the cloud's cross-section and 9 the polar angle, 
measured from the normal to the radiating disc. The 
anisotropy arises from the change in projected surface 
area of the infinitesimally thin accretion disc of bolometric 
luminosity L assumed here. Thus, the cloud is moving in a 
force field (per unit of mass) 



f( a\ GM 

/M) = — 



31 



cos(0)| - 1 



(4) 



where fi is the mean molecular weight and N c i — \mj ' ^m v a 
is the column density through the centre of a single spherical 
cloud, which is equivalent to an optical depth and assumed 
to be constant throughout the evolution. I = L/L c dd is 
the Eddington ratio where L 0( jd = 4ttGM m p c/aT is the 
Eddington luminosity and ot the Thomson cross-section. 
It is straightforward to see that - in general - the radiation 
field effectively changes the gravitational constant G or the 
central mass M continuously throughout the orbit. Angular 
momentum conservation keeps the orbital plane fixed in 
space and hence equatorial orbits are unaffected by radiation 
pressure. We discuss possible contributions from other forces 
like, e. g., drag forces in Sect. [6] 

3 CLOUD ORBITS 

Since the total force on a cloud is always radial, specific 
angular momentum L is conserved and thus the orbit 
remains in a plane with fixed inclination i against the 
accretion disc. It is sufficient to consider ^ i ^ 7r/2. Let ip 
be the angle in this plane as measured from the ascending 
node, then cos(0) = sin(i) sin(0) and the radial equation of 
motion is 

.. L 2 



GM 



(k\ sin(^)| - 1 



where L = r 2 ip and 
■ sin(j). 



(5) 



(6) 



LlN c l(TT 

Changing the independent variable from time to angle ip and 
substituting u = 1/r, the equation of motion can be written 

as 

d 2 u 



_ GM 

diP 2+U ~ L 2 ' 



k\ sin(^)| 



This differential equation has the solution 



u(ip) 



GM 



k 

1 + — ip cos(i/>) + a sm(i/>) + /3 cos^) 



for ^ if) ^ 7r and similarly for ir < ip ^ 27r, 



GM 



1 — — ip cos(-0) + (a + k) sin(?/') + 



(7) 



(8) 



(9) 



+ (13 + kn) cos(-i/>) 



with arbitrary constants a, /3 and under the condition that 
both u and du/dip are continuous at ip — n. The parameters 
M and L appear as scaling factors and have no influence on 
neither stability nor shape of the orbit. 



4 STABILITY OF ORBITS 

The orbit equations [S] and allow for unbound orbits 
in such way that u can have a root on ^ ip ^ 2n. 
Requesting it to be at least positive at ip = n ■ n/2 for 
all integer n, results in the two constrains — 1<q<1 — k 
and —1 < P < 1 — kir/2. The latter hints that there is a 
critical value of k c — 4/tv, for which bound orbits become 
impossible. This has been confirmed numerically by solving 
the respective orbits for an extensive number of parameter 
combinations. Thus, equation [6] implies a critical inclination 
for any particular choice of column density and luminosity, 
given by 

' 4 fj,N c i<jT ' 



i c = arcsm 



V TV 



31 



(10) 



above which all orbits are unstable. 

Some typical values are illustrated in Fig. [T] where 
H — 0.61 is chosen representative for solar metallicity gas. 
For example, in a typical Seyfert galaxy with Eddington 
ratio of 0.1, a cloud with column density 10 23 cm~ 2 is 
restricted to angles smaller than 10° . If the column density is 
greater by a factor of about 6, polar orbits become allowed. 
A similar argument can be made if the column density is 
kept fixed, namely that the greater the Eddington ratio of 
the source, the smaller the region where stable orbits are 
possible. 



90° 
80° 
70° 
60° 
50° 
40° 
30° 
20° 
10° 



N d = 10 22 cnr ; 

N d = 10 23 cm- : 

N d = 10 24 cm- ; 



0. 0.1 0.2 0.3 0.4 0.5 0.6 
I 



0.7 



0.9 1. 




Figure 1. Critical orbit inclination vs. Eddington ratio of the 
central source and cloud column density (fi = 0.61). 
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Figure 2. Distance of the cloud to the central mass over one 
orbit for case (e) from Fig. \3\ 



Because k c > 1, the critical inclination is always larger than 
the inclination defined by force equilibrium (/ = 0), up to a 
factor of ~ 1.74. This is plausible, because clouds on orbits 
with a comparatively high inclination spend only a relatively 
short amount of time in the radiation-dominated region. 

5 SHAPE OF STABLE ORBITS 

It can be verified that all bound orbits are closed by 
evaluating u and du/dip at ip = e from equation [H] and at 
ip = 2-7T — e from equation [9] and letting e go to zero. That is 
to say the total energy averaged over one orbit is constant, 
although the force field is non-conservative. Closed orbits are 
not expected under different assumptions, such as variable 
N c i or I, which might be expected for a pressure-confined 
cloud with constant mass adjusting instantly to pressure 
equilibrium within a stratified atmosphere. 

The family of stable orbits features a range of different 
orbit shapes. For small values of k, orbits are close to 
elliptical and reduce to Kepler orbits if k approaches zero. 
This is the case for either zero inclination, zero luminosity 
or very high column density. If the value of k is close to 
critical, all stable orbits are eccentric, show a waist and can 
also be asymmetric. There are then two local maxima of 
the cloud's distance to the central mass over one orbit, as 
shown for example in Fig. [2] In addition to the eccentricity 
e calculated from apocentre and pericentre of the orbit, it is 
convenient to use another similarly defined quantity 



r(V>i) - r(jfe) 
r(V>i) + r{ip2) 



(11) 



calculated from these two maxima, together with the 
angular offset Aip = ip2 — tpi between them, to distinguish 
between characteristic shapes. 

Some exemplary orbits are shown in Fig. [3] with 
parameters summarized in Table [T] The transition from 
orbit (a) to (e), which have identical k, can be understood in 
such way that the near elliptical orbit (a) is the most circular 
orbit possible in this case. Orbit (b) has two distinguishable 
apocentres (e' 7^ 0) and shows a pronounced waist, but the 
apocentres still lie opposite to each other {Aip = n). In 
contrast to that orbit (c) is completely asymmetric, which 
is the most prevalent shape. There is again an axis of 
symmetry for orbit (d) because here e' vanishes, although 







(b) 





Figure 3. Characteristic shapes of orbits depicted in the orbital 
plane according to equations \E\ and [9] Distances are marked in 
units of GM/L 2 . The parameters for each orbit can be found in 
Tabled] 



Table 1. Parameters for the example orbits in Fig. [3] 
illustrated in section [5] 





k 


a 


P 


e 


e' 


Aip - 


- 7T 


(a) 


1 


-0.5 


-0.785 


0.40 










(b) 


1 


-0.5 


-0.65 


0.73 


0.63 







(c) 


1 


-0.25 


-0.65 


0.90 


0.75 


0.32 


7T 


(d) 


1 


-0.25 


-0.785 


0.63 





0.26 


7T 


(e) 


1.2 


-0.4 


-0.935 


0.91 


0.21 


0.16 


7T 



Aip 7^ 7T. Orbit (e) is of similar type as orbit (c), but for 
a slightly higher k, which emphasizes the above-mentioned 
waist . 

6 DISCUSSION 

This Letter represents a first step towards understanding 
the dynamical properties of optically thick clouds in point 
potentials including an anisotropic radiation source. Such 
simulations and theoretical predictions can be used to assess 
the dynamical structure as well as the stability of AGN tori 
and the BLR. 

For a general treatment of the problem, several 
simplifications are necessary. First of all, we assume that 
the clouds move through vacuum. In reality, e. g., AGN tori 
are thought to be made up of a multiphase medium, where 
dense and cold clouds are embedded in a hot and diffuse 
ionized gas component. This intercloud medium leads to ram 
pressure interaction with the cloud, which has two effects: 

(i) locally it compresses the front part of the cloud and 

(ii) globally it leads to a deceleration and thereby a change 
of the orbital properties and finally an inspiral of the cloud 
towards the centre. Related is the is sue of confinement of th e 
clouds, as hydro instabil ities (e. g. ICowie fc McKeelll977r ). 
radiation pressure (e. g. ISch artmann ct al. Il201lh. thermal 
conductivity (e. g. iMcKee fc Ostrikerl Il977l ; iBurkert et al.l 
l2012l;ISchartmann et al.ll2 012) or collisions with other clouds 
IjKrolik fc Begelmanlll988l ) would lead to the dissolution of 
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the clouds on time-scales of the order of the dynamical 
time or even faster, depending on the geometry of the 
cloud distribution. An additional issue is the built-up of a 
dense underlying disc through dissipative collisions, which 
would eventually lead to a collapse of the cloud distribution. 
Such equatorial discs are indeed found in coexistence 
with geometrically thick AGN tori in nearby Seyfert 
nuclei, probed with the help of near-i nfrared interfer o metry 
as well as in maser emission ([Tristram et al.l 120071 : 
iRaban et ail 120091 ; ISchartmann et alj |2010| ). Sufficiently 
elastic cloud-cloud collisions (e. g. due to strong internal 
magnetic fields) would enable the exchange of energy and 
thereby alter the initial orbital structure as well, leading 
at the same time to an accretion flow towards the centre 
and an in crease of the veloc i ty di s persion of the clou d 
ensemble (|Krolik fe Begelmanl 1 19881 ; Ivollmer et ail |2004) . 
Depending on the confinement, clouds might change their 
size when, e. g., moving through a stratified atmosphere 
which then in turn affects the ra diation pressure force and 
hence the orbital e volution (e. g. iNetzer fe Marziamll201ol ; 
iKrause et al.l 1201 ll ) . In addition magnetic fields influence 
the dynamical evolution and might contribute either to a 
faster destruction o r to long lasting confinement of clouds 
l|Krause et al.ll2012T ). depending on the - generally unknown 
- magnetic field topology. 

7 CONCLUSIONS 

The shape and stability of an optically thick cloud's orbit 
around a central black hole with a surrounding accretion 
disc are strongly affected by radiation pressure. Although 
deviating from ellipses, we show that all bound orbits are 
closed. The stability of an orbit is determined by the ratio of 
the disc's Eddington fraction to the cloud's column density 
and in particular the inclination of the orbit. 

This has a possible application in the study of AGN 
tori, which are composed of dusty clouds of gas. The clouds 
can reasonably be assumed to have a high optical thickness 
because dust absorption and the emitted spectral energy 
distribution of the central engine both peak in the UV 
regime. Therefore, radiation pressure is likely to influence 
the dynamics of clumpy AGN tori significantly. Compared to 
the derived limit on the inclination of stable cloud orbits, the 
geometrical thickness of a real torus might be amplified due 
to an outflow consisting of ejected clouds (lEmmering et al.l 
Il992l ; iKonigl fe Kartidll994l ; lElitzur fe Shlosmanll2006l . and 
references therein) . These could originate from a distribution 
with large covering angle at the onset of activity or from 
an accretion flow triggered by cloud-cloud collisions. This 
outflow would be of particular importance for large ratios 
l/N c i, e. g. in Quasars. When supply for the outflow has 
ceased, there could still be left a rather thin but long-lived 
torus. 

However, for a full treatment of cloud ensembles, the 
interaction of multiple clouds would need to be considered, 
such as fluctuating radiation pressure caused by shadowing 
of the central light source, secondary radiation pressure 
from heated dust and merging or collisions of clouds. In 
this environment, magnetic fields might be able to support 
clouds against destruction by tidal or hydrodynamic forces 
to ensure the stability of clouds on the time-scale of one 
orbit and also provide the necessary coupling of the dust 
and gas components. 
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